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' The rotating neutron superfluid in the inner crust of a neutron star is threaded 

by quantized vortex lines. The pinning force from lattice nuclei and the Magnus 
force from neutron superfluid act onto a vortex line that has a finite tension. 

\q • The configuration of a vortex line in equilibrium comprises a number of pinned 

straight lengths separated by unpinned kinks when the rotation axis is slightly 

Q\ • inclined to the major axis of a crystal lattice and the Magnus force is not so strong. 

^ 1 Energy of ~ 4 MeV is required to form a kink at densities of ~ 3.4 x 10 13 g cm -3 . 

Magnus force makes the vortex line parabolic on global scale. There exist two 

i-C ■ modes, rotational and helical (Kelvin), for oscillations allowed on a vortex line. 

The vortex oscillations in the pinned straight segments are possible only above 
a minimum frequency. We find no unstable mode that grows with time. Hence, 
the vortex configurations with kinks may be stable and yield the strongly pinned 
state. The essence of our results may not be altered even when we consider the 

J^j ■ vortex behaviors in a polycrystalline structure. Our studies suggest that pinning 

of a vortex line in a polycrystalline structure would be still strong enough to 
explain the large glitches. 

Subject headings: dense matter — pulsars: general — stars: interiors — stars: 
neutron 



1. INTRODUCTION 

Glitches, observed from more than 30 isolated pulsars (Lyne, Shemar, & Smith 2000; 
Wang et al. 2000), are a sudden jump in pulsar rotation rate. Slow recovery of a rotation rate 
following these events implies the presence of loosely-coupled fluids or neutron superffuids in 
the interior of a neutron star. The transfer of angular momentum from the rotating neutron 
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superfluid to the normal matter is likely to be responsible for pulsar glitches and postglitch 
relaxation. 

A rotating neutron superfluid is threaded by quantized vortex lines. The superfluid 
can alter its angular velocity by the radial motion of vortices. In the inner crust of a 
neutron star neutron rich nuclei that form the bcc lattice coexist with neutron superfluid. 
The vortex lines in the inner crust are normally pinned to the lattice nuclei. Difference in 
velocity between superfluid and nuclear lattice builds up since magnetic braking slows down 
the nuclear lattice and pinning prevents the vortex lines from moving. Then, the neutron 
star undergoes sudden unpinning of a large number of vortex lines followed by the outward 
motion. This catastrophic unpinning has long been considered as a promising cause for 
glitches (Anderson & Itoh 1975; Ruderman 1976; Alpar et al. 1984; Pines & Alpar 1985; 
Epstein, Link, & Baym. 1992; Baym, Epstein, & Link 1992). 

A vortex line is subjected to the Magnus force when it moves relative to the superfluid. 
The catastrophic unpinning model assumes that the pinning force is strong enough to sustain 
the vortices to the nuclear lattice against the Magnus force until the moment just before 
glitches. Based on the condensational and kinetic energies, the pinning energy is estimated, 
ranging approximately 1 — 10 MeV per nucleus (Alpar et al. 1984; Epstein & Baym. 1988; 
Pizzochero, Viverit, & Broglia 1997). This magnitude itself is strong enough against the 
Magnus force expected just before glitches. Jones (1992) mentions, however, that the pinning 
forces on a randomly oriented rigid vortex line largely cancel since there are nearly an equal 
number of pinning sites on either side of the vortex line. Link, Epstein, & Baym (1993) 
find that for a vortex line of finite tension, pinning becomes much more efficient by slightly 
bending and forming kinks. Recently, Jones (1997, 1998, 1999) argues that vortex interaction 
with a polycrystalline structure does not provide pinning strong enough to explain the large 
glitches observed in the Vela pulsar (see Section 5 for detail). 

Here we study the vortex configurations, oscillations and pinning in the inner crust of a 
neutron star. In Section 2 we derive the equation of motion of vortex lines. In Section 3 we 
present the equilibrium configurations of vortex lines. In Section 4 we examine the stability 
of equilibrium configurations and the oscillations excited on a vortex line. In Section 5 we 
discuss the vortex pinning based on the results of the previous sections. In the last section 
we summarize the results and mention the conclusions. 
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2. EQUATIONS OF MOTION OF VORTEX LINES IN AN INNER CRUST 

We consider the motion of vortex lines in an inner crust, where the lattice nuclei are 
embedded in a superfluid neutron sea. We adopt the Cartesian coordinates fixed in the 
rest frame of the lattice nuclei. The z-axis is chosen to be parallel to one of the major 
axes of a body-centered crystal lattice. In this paper the most attention will be devoted 
to configurations and oscillations of vortex lines which are parallel to or slightly inclined 
against the z-axis. We limit the inclination angle of a vortex line against the z-axis to a 
small value. Then, we can describe their motion and configurations by the two-dimensional 
position vectors, u = (x, y) , of vortex lines in the xy plane which depend on the coordinate 
z. The equation of motion of vortex lines is written as 

f du \ ^d 2 u , , 

f r + pKX (-- v ,)+T^ = , (!) 

where p is the superfluid density, f p the pinning force per unit length of a vortex line, k the 
circulation of a vortex, v s the locally uniform superfluid velocity and T the tension of the 
vortex line. 

Vortices in the superfluid are quantized. Each vortex line carries one quantum k = h/2m 
of circulation, where h is the Planck constant and m the neutron mass. Since we consider 
the vortex lines aligned or almost aligned with the z-axis, we approximate the circulation 
vector by 

h „ 

where z is the unit vector along the z-axis. 

In equation (1) the second term represents the Magnus force which arises when there 
exists the relative motion between the local superfluid and the vortex. The magnus force 
acting on a unit length of a vortex line is given by 

f M = pax (v v - v a ) (3) 

where v v is the vortex velocity. In the vortex unpinning model the difference in angular 
velocity between the inner crust superfluid and the crust is reduced significantly just after 
the glitches, decoupling the inner crust superfluid from the crust. As the crust is slowed 
down by the external braking torque, the angular velocity difference Afl builds up between 
the local superfluid and the lattice nuclei. Using the glitch interval t g and the angular 
deceleration rate Cl, the angular velocity difference just before the glitch is estimated as 
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AQ ~ \n\t g . (4) 

In the Vela pulsar the angular velocity difference and the Magnus force acting on the pinned 
vortex line may amount to ~ 0.01 rad s^ 1 and ~ 6.8 x 10 14 dyne cm -1 , respectively, just 
before the glitches. 

The third term in equation (1) comes out from the tension of vortices when vortices are 
distorted from the straight lines. The energy of vortex lines is equal to the kinetic energy of 
the velocity field induced by them. The tension is given by the energy per unit length, 



T = pnc K , (5) 

where the parameter ck is only logarithmically dependent on the wave number and the 
neutron superfluid coherence length (Sonin 1987). For numerical estimates we assume 
c K = 1.2 x 10~ 3 cm 2 s- 1 (Jones 1998). 

The vortex-nucleus interaction has an attractive component arising from the change in 
the condensation energy and a repulsive component arising from the change in kinetic energy 
of circulating superfluid. We adopt only the attractive component since we are interested in 
the inner portion of the inner crust, where the attractive component is dominant over the 
repulsive component. The interaction potential we adopt is 

V = -E p ew(-^, (6) 

where E p is the pinning energy (or the interaction energy) per nucleus, £ is the vortex core 
radius or the coherence length of the neutron superfluid and u is the distance in the xy plane 
between the nucleus and the vortex line, u 2 = x 2 + y 2 . The coherence length of the neutron 
superfluid is written as 

£ = — ■ (7) 

where A is the energy gap for the superfluid neutrons, Ep their Fermi energy and hp the 
Fermi surface wave number (Alpar 1977; Alpar et al. 1984). The pinning force, which acts 
to pin the vortex line to nuclei, per site is approximately given by 

F P = ^. (8) 



We can obtain the lattice-vortex interaction potential by summing the above potential 
(6) contributed from each nucleus over the entire nuclei of a single large crystal. Note that 
the lattice-vortex interaction potential has the translational properties and hence can be 
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expanded into the Fourier series using the reciprocal lattice vectors. We assume the lattice- 
vortex interaction potential to be uniform in the z-direction since we are interested in the 
configurations and oscillations of vortices, which are approximately aligned with the z-axis, 
with the scale and wave lengths significantly larger than the lattice constant. Hence, of the 
twelve basic reciprocal lattice vectors of the bcc lattice we adopt the basic four vectors in 
the xy plane, (g, g), (—g, g), (—g, —g) and (g, —g), where g = 2-n/a with a being the lattice 
constant. Then, the lattice- vortex interaction potential per unit length of a vortex line is 
given by 

U p = 2U g [cos{gx + gy) + cos{-gx + gy)}, (9) 

where the Fourier coefficient U g is ~ — 9.7 x 10 3 erg cm -1 at densities of p ~ 3.4 x 10 13 g cm~ 3 . 
We obtain the pinning force acting on the unit length of a vortex line by taking the derivative 
of Up with respect to x and y, 

f px = 2gU g [sm(gx + gy) - sin(-gx + gy)] (10) 

and 

f py = 2gU g [sm(gx + gy) + sin(-gx + gy)}, (11) 

where f px and f py are the x and y components of the pinning force, respectively. 

We conduct the coordinate rotation and transform the variables from (x, y) to (0, tp) 
defined as 



(f) = gx + gy (12) 

and 

^ = -gx + gy- (13) 
Using equations (10), (11), (12) and (13), the equation of motion of vortices (1) reduces to 



2 TT ■ th -' , , i)2 ' 



and 



Ag U g sin <p - pn— + pngv s + pncx—^ = (14) 
at oz z 



2 d4> d 2 ip 

Ag UgSmip + pK— - pngv s + pnc K ^r -j = 0. (15) 



dt ' » b ' r n dz 2 

Here and hereafter we assume for simplicity that the neutron superfluid flows in the y- 
direction. Note that equations (14) and (15) describe the motion of vortices in terms of the 
components in the directions of reciprocal vectors. 
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3. EQUILIBRIUM CONFIGURATIONS OF VORTEX LINES 



3.1. Equations for Equilibrium Configurations 



If the vortex lines are completely rigid, the vortex configuration would be simply a 
straight line. The tension of a superfluid vortex line in the inner crust of a neutron star, 
however, is finite. Hence, when the vortex lines undergo the forces such as the pinning and 
Magnus forces, the vortex lines in general are bent and their shape is deformed from the 
straight line. We examine the equilibrium configurations of vortex lines here. 

We assume that in equilibrium dtfi/dt = and dip/dt = in equations (14) and (15) 
and hence the vortices are stationary with respect to the nuclear lattice. Then, equations 
(14) and (15) can be integrated to yield 



Once a displacement and orientation of the vortex line are given at , for example, z = 0, we 
can solve equations (16) and (17) and derive the equilibrium configurations of vortex lines. 

Hereafter, for numerical calculations we use the physical parameters such as the pinning 
energy, the coherence length and the lattice constant at densities of p ~ 3.4 x 10 13 g cm~ 3 , 
where a large fraction of crust superfluids resides and the pinning energy peaks. We adopt 
v s ~ 10 4 cm s -1 for the superfluid velocity unless otherwise stated. 

The first, second and third terms in equations (16) and (17) correspond to the pinning 
force, the Magnus force and the vortex tension, respectively. The pinning force per site 
is F p ~ 9.6 x 10 5 dyne and the vortex tension is T ~ 8.1 x 10 7 dyne, while the Magnus 
force acting on the vortex section of a length corresponding to the lattice constant a is 
pnv s a ~ 3.9 x 10 3 dyne. We note that the Magnus term is much smaller than the pinning 
and tension terms. 

We should note that equations (16) and (17) are analogous to the energy conservation 
equation for the particle motion in a potential well. The first and second terms and third 
term on the left hand side of equations (16) and (17) may correspond to the potential and 
kinetic energies in the particle motion, respectively, while the integration constant on the 
right hand side to the total energy. The gradient of the vortex line against the z-axis may be 




(16) 



and 




const. 



(17) 
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looked as the velocity of a particle. The difference between the "total and potential energies" 
yields the gradient of the vortex line. The pinning term in equation (16) (equation (17)) 
has peaks at = 1m\ (ip = 2mr) with n being integer, while the Magnus term increases 
(decreases) linearly with (ip). As easily inferred from this analogy, solutions to equations 
(16) and (17) can be classified into the three characteristic cases in terms of the magnitude 
of the integration constant ( the "total energy") relative to the peak values of the potential 
energy due to the pinning and Magnus forces. 



3.2. Straight Line Solutions 

Let us consider the case where the integration constant (the "total energy") in equations 
(16) and (17) is just equal to a peak value of the sum of pinning and Magnus terms (the 
"potential energy"). The derivative of <f> and ip with respect to z should be zero over an 
entire vortex line. Hence, the equilibrium solution for the vortex configuration is a straight 
line which is in parallel with the z-axis or the major axis of a single crystal lattice (Fig. 
1). Note that the average orientation of a vortex line expresses the direction of the rotation 
axis. Thus, the rotation axis of the star is also parallel to the z-axis. When the Magnus 
term is neglected in equations (16) and (17), the straight vortex line passes through the 
bottom of the pinning potential well along the z-axis. When the Magnus term is included, 
the straight vortex line, though still parallel to the z-axis, is displaced only slightly from the 
z-axis, x ~ 1.4 x 10 -14 cm and y ~ cm. In actual cases when the Magnus force acts on 
the vortex lines, however, the equilibrium configuration of the vortex may deviate from the 
straight line, though very slightly, since the pinning sites are located discretely along the 
z-axis. We have the strongest pinning regime when the rotation axis happens to be oriented 
along the major axis of a single crystal since the vortex line intercepts the pinning sites most 
often. 



3.3. Sinusoid-like Solutions 

Let us move to the case where the integration constant (the "total energy") in equations 
(16) and (17) is smaller than peak values of the sum of pinning and Magnus terms (the 
"potential energy"). This case corresponds to the motion of a particle which is bound in the 
potential well. As the particle exhibits a periodic motion, the vortex line goes back and forth 
between the two neighboring lattice planes (Fig. 2). The average direction of the vortex line 
and hence the rotation axis of the star are parallel to the z-axis. Note that the significant 
energy is required to deform the vortex from the straight line to the sinusoid-like shape ( see 
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the next sub-section for the calculation of the total energy of a vortex). Hence, it may be 
the straight vortex line which is realized in nature when the rotation axis is parallel to the 
z-axis or the major axis of a single crystal lattice. 



3.4. Kink Solutions 

For simplicity and clarity we neglect the Magnus force here. Let us consider the case 
where the integration constant (the "total energy") in equations (16) and (17) is larger than 
the peak value of the pinning term (the "potential energy"). This case corresponds to the 
motion of a particle which continues to move in one direction with varying speed. The 
derivative of and tp with respect to z is not now zero everywhere along a vortex line. 
Hence, we examine here the equilibrium vortex configuration when the average orientation 
of a vortex line and hence the rotation axis of the star are inclined against the z-axis (or the 
major axis of a single crystal). 

We study first the equilibrium solutions when the angle between the rotation axis of a 
star and the major axis of a crystal is infinitesimally small. Over an entire interval of a large 
single crystal the vortex passes from one lattice plane to the adjacent lattice plane. Hence, 
we set the boundary conditions such as = at z = — oo and = 2n at z = oo. Then, we 
can solve equations (16) analytically and derive the equilibrium configurations, 



= 4 arctan 

where 



cxp 



(t)] ■ <«■> 



Figure 3 shows that the vortex can move from one lattice plane to the adjacent lattice plane 
by bending itself and forming a kink when the vortex tension is finite. The vortex line is 
displaced substantially from the pinning sites in the kink part, while in the other part the 
vortex line is straight and passes through nearly all of the pinning sites along the z-axis. I 
given by equation (19) expresses a length of the kink. In the middle of a kink the vortex line 
is inclined against the z-axis by 2.5°. 

Energy is needed to create a kink from a straight vortex line since the length of a vortex 
line increases and a part of it should overcome the pinning potential. The self-energy of 
vortex lines is equal to the kinetic energy of the velocity field induced by them. The increase 
of the vortex self-energy due to bending is 
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dz, 



(20) 
(21) 



while the energy necessary to overcome the pinning potential is 



AE P = / 2U g (cos(f> + cos ip)dz- 2[/ 9 (cos0 + cos0)dz. (22) 

J —oo J —oo 

The length of a kink can be determined from minimizing the sum of AE S and AE P (Link 
& Epstein 1991; Ruderman 2000). Using the solution (18) together with if) = 0, the total 
energy necessary to form a kink is calculated as 

AE t = AE S + AE P = 8 ^> KC * ) 2 , (23) 

which yields ~ 4.0 MeV at densities of ~ 3.4 x 10 13 g cm -3 . 

We increase gradually the integration constant in equations (16) and (17) and hence 
the inclination angle of the rotation axis against the major axis (z-axis) of a single crystal. 
The vortex lines now should thread through many lattice planes. We find that when the 
inclination angle is small and less than 0.7°, vortex lines in equilibrium configurations move 
from one lattice plane to the next by forming a clear kink and the number of kinks increases 
with increasing inclination angle (Fig. 4). The separation L z in z between the neighboring 
kinks varies as L z oc 6~ l . The property of a kink in Fig. 4 is the same as that in Fig. 3. A 
vortex line can orient its average direction to the direction of the rotation axis by adjusting 
a number of kinks. Apart from the kink parts, the vortex lines are straight and in parallel 
to the z-axis, and are firmly pinned to the lattice nuclei. Even though the extra energy is 
required to create kinks, the total energy of a vortex line with kinks is much lower than that 
of the inclined straight vortex line without kinks since except for kink parts the vortex line 
lies in the bottom region of the pinning potential (Link, Epstein, & Baym 1993). 

As the average orientation of a vortex line inclines against the major axis, the separation 
between the neighboring kinks and hence the portion of a vortex line available for pinning 
decrease. As seen in Fig. 4, the straight pinned part of the vortex line becomes shorter than 
the kink part when the inclination angle is above ~ 0.7°. The vertical pinned portion as well 
as the kink feature becomes less prominent at ~ 1.3°. When the inclination angle exceeds 
~ 2.5°, the equilibrium configuration of a vortex line becomes almost the straight line. This 
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consequence arises from the fact that the portion of a vortex line for pinning becomes too 
short to compensate the energy cost for producing kinks when the inclination angle exceeds 
a certain value. 

When the average orientation of a vortex line inclines significantly against the z-axis, 
we need to include the z-components for both the force and displacement, though neglected 
here. We expect that the lattice planes different from those considered here may become 
important for pinning. The three dimensional calculations are needed in order to know the 
equilibrium configuration of a vortex line accurately. 



We now examine how the kink solutions are modified when we include the Magnus 
terms in equations (16) and (17). The pinning and Magnus terms as a function of <fi (ip) 
are represented by the cosine curve and the inclined straight line, respectively. Hence, there 
exists 0c (ip c ) at which the integration constant (the "total energy") becomes equal to the 
sum of the pinning and Magnus terms (the "potential energy"). The integration constant is 
larger (smaller) than the sum of the pinning and Magnus terms at < <p c (ip > ip c ), while 
smaller (larger) at <fi > <j) c (ip < ip c ). This case corresponds to the motion of a particle which 
moves forwards and then at the potential wall bounces to move backwards. 

Figure 5 illustrates the equilibrium configurations of vortex lines for different superfluid 
velocities v s relative to the lattice nuclei. The equilibrium configuration is characterized 
locally by a kink and globally by a parabola. As explained above, the kink arises from the 
balance between the pinning and tension terms. As expected from the analogy with the 
particle motion, the vortex lines extends first rightwards and then leftwards as z increases. 
When we take a long scale average of equation (1), the pinning forces, which alternate 
the sign (the direction) along a vortex line, cancel out whereas the Magnus force, which is 
constant along a vortex line, is left. Hence, when combined with the tension term in equation 
(1), the constant Magnus term yields the parabolic configuration. The curvature radius of a 
vortex line is given approximately by 



As the superfluid velocity relative to the lattice nuclei increases, the curvature and hence 
the deformation due to the Magnus force becomes larger (Fig. 5). Notice also that the kink 
feature becomes less prominent as the deformation becomes larger due to the same reason as 
mentioned in the previous sub-section. The three dimensional calculations are needed again. 



3.5. Effect of the Magnus Force 




(24) 



3.6. Effect of a Density Variation 



Numerical estimates in this paper are made mostly by adopting the physical parameters 
at the inner part of an inner crust, where the density is higher and the pinning energy peaks. 
The kink properties and pinning strength depend on matter densities through the nature of 
nuclei and superfluids. For example, the kink length and creation energy are l/a ~ 7.3 and 
AE t ~ 4.0 MeV at densities of p ~ 3.4 x 10 13 g cm -3 in the inner part of an inner crust, 
while l/a ~ 39 and AE t ~ 0.089 MeV at densities of p ~ 2.7 x 10 12 g cm" 3 in the outer 
part of an inner crust. The choice of parameters may affect the results. Note, however, that 
most of crust superfluids reside in the inner part of an inner crust. Hence, even if the outer 
part of the inner crust is included adequately, the change may be small and may not alter 
the essence of our conclusions. 



3.7. Effect of a Polycrystalline Structure 

The inner crust nuclei are likely to take a polycrystalline structure (Jones 1998). The 
orientation of a crystal lattice may change from one domain to another. Then, along a 
vortex line there may lie a large number of single crystals with various orientations against 
the average direction of a vortex line. The important physical parameters which vary with 
the orientation of a crystal lattice are the distance between the neighboring lattice planes 
for pinning and the distance between the neighboring pinning sites along a vortex line. It is 
likely that even when the vortex line is not close to the main axis of a crystal lattice, the 
equilibrium configuration of a vortex line may also be expressed by a kink solution with the 
pinning lattice planes different from those containing the main axis of a crystal lattice. As 
seen from equations (19) and (23), the kink length and creation energy are in proportion 
to the distance between the neighboring lattice planes for pinning and depend weakly on 
the Fourier coefficient (U g ) of a pinning potential. The kink properties may vary, though 
not significantly, depending on the orientation of a crystal lattice. The pinning strength 
of a vortex line is inversely proportional to the separation of the neighboring pinning sites 
along a vortex line. When the vortex line is not close to the major axis of a crystal lattice, 
the pinning site distance will be longer and hence the pinning of a vortex line will be less 
strong. The various orientations of crystal lattice along a vortex line may reduce the pinning 
strength of a vortex line as a whole compared to the case with a small inclination angle. 
We expect, however, that pinning, though weaken, is still strong enough to explain the large 
glitches since the change of the pinning site distance more than an order of magnitude is less 
likely. We need to conduct the three dimensional calculations in order to understand the 
equilibrium configuration and pinning strength accurately for various orientations of crystal 



lattice as well as for large Magnus forces. 



4. VORTEX OSCILLATIONS 



Let us now consider the oscillations excited on and propagate along the vortex line 
(Link, Epstein, & Baym 1993). In order to examine oscillations of small amplitude about 
equilibrium, we write 

= 0o + 0i (25) 

and 

V> = Vo + V>i, (26) 

where the subscripts and 1 denote the equilibrium solutions and the perturbations, respec- 
tively. Using equations (25) and (26), we linearize equations (14) and (15) and obtain 

Ag 2 U g cos 0o x 0x - p K -±± + Pkck-^- = (27) 

and 

Ag 2 U g cos x^i + ptS^ + pnc K ^- = (28) 

at az 1 

We consider plane- wave solutions oc exp(ipz — iut) for 0! and ipi, where uj is the angular 
frequency and p is the wave number. Then, we derive the dispersion relation for oscillations, 

(Ag 2 U g cos 0o — pkckV 2 ) (^9 2 U g cos -0o — PkckP 2 ) — (pnuj) 2 = 0. (29) 

If we adopt such solutions as ipo = — 0o for equilibrium, from equation (29) we obtain the 
dispersion relation simplified as 



2 , y c/ } 

crP H — COS 0o 

pK 



(30) 



Next, we examine the property of oscillations represented by the dispersion relation 
(30). The first term in the dispersion relation (30) can be derived when in equations (27) 
and (28) we include only the second and third terms ( the Magnus and tension terms) 
ignoring the first term (the pinning term). Hence, the first term in the dispersion relation 
(30) expresses the Kelvin wave, which is circularly polarized and propagates helically along 
a vortex line (Sonin 1987). The second term in the dispersion relation (30) can be derived 
when in equations (27) and (28) we include only the first and second terms (the pinning and 
Magnus terms) neglecting the third term (the tension term). We find that the second term 
in the dispersion relation expresses the circular motion of a vortex line as a whole around 



the equilibrium position. In general the dispersion relation (30) represents the combined 
oscillations of Kelvin and rotational modes. 

Figure 6 illustrates the dispersion relation given by equation (30) for different O - The 
angular frequency of allowed oscillations is plotted as a function of wave number. At smaller 
wave numbers the angular frequency of oscillations is determined mainly by the rotational 
mode, while at larger wave numbers mainly by the Kelvin mode. We should note that there 
exists the minimum frequency for the oscillations especially when O is close to 2nn. The 
minimum frequency for 0o — 2wr is 

u - 

pn 

which yields 7 x 10 17 rad s^ 1 at densities of ~ 3.4 x 10 13 g cm -3 and 1 x 10 16 rad s' 1 at 
densities of ~ 2.7 x 10 12 g cm -3 . O ~ 2nir in the pinned segments, where the vortex line lies 
deep in the pinning potential well. On the other hand, at the kink segments where 0o ~ (2n+ 
l)ir, the minimum frequency is zero. These results lead to the important conclusion that the 
vortex oscillations with frequencies lower than the minimum frequency, uo < 7 x 10 17 rad s _1 , 
can not be excited on and can not propagate in the pinned part of a vortex line, whereas 
oscillations with frequencies higher than the minimum can propagate along a vortex line, 
varying the wave length. The vortex oscillations with lower frequencies are allowed in the 
kink part of a vortex line. We should note, however, that the dispersion relation (30) in the 
kink part is less reliable for lower frequencies since the wave length exceeds the kink length 
and hence the local perturbation analysis is not adequate. We need to conduct the non-local 
analyses for the accurate dispersion relation in low frequencies. 

In a polycrystalline structure the vortex line passes through a number of single crystals 
with various orientations. The local perturbation analysis can be applied to each single 
crystal to examine the vortex oscillations. Equations (30) and (31) may also be derived 
for each single crystal, although the physical parameters in one domain differ from those in 
another. We can expect the physical parameters to change rather sharply at the boundaries 
between the single crystals. If the frequencies are sufficiently high above the minimum 
frequencies, the oscillations will propagate along a vortex line passing through many single 
crystals. If the frequencies are around the minimum frequencies, the oscillations are reflected 
at the boundaries and in some cases will be trapped within one crystal or some. If the 
frequencies are lower than the minimum frequencies, there may be no oscillations which 
propagate over a number of single crystals. 



(31) 
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5. VORTEX PINNING 

It has been considered that pulsar glitches are a consequence of the large scale unpin- 
ning of superfluid neutron vortices from the lattice nuclei in the inner crust of neutron stars. 
Recently, however, Jones (1997, 1998, 1999) argues that vortex interaction with a poly- 
crystalline structure does not provide pinning strong enough to explain the large glitches 
observed in the Vela pulsar. Let us examine this conclusion in the light of our new results 
mentioned in the previous sections. 

Jones describes the reason for the conclusion that the structural regularity of a defect- 
free single crystal allows displacement of a vortex, under a Magnus force, to a continuous 
sequence of new positions with no more than extremely small changes in energy. Further- 
more, Jones considers the motion of kinks propagating along a vortex line with a velocity of 
Vk- For simplicity we consider the vortex lines which have kinks bent in the y-direction and 
extend approximately in the z-direction. The kink propagation along a vortex line yields a 
displacement and hence motion of the vortex line in the y-direction. If the kink propagation 
velocity is 

v K = —v s , (32) 
a 

the time average of a vortex velocity in the y-direction is equal to the superfluid velocity 
(Jones 1997, 1998). In effect, vortex lines move in the y-direction together with the superfluid, 
almost not suffering from pinning. These energy and velocity arguments lead Jones to the 
above conclusion that pinning would be extremely weak and insufficient for the large glitches. 

We point out several difficulties included in these arguments. Kinks should be created 
and supplied around the edge of a single crystal at a rate of one kink every a/v s seconds 
in order for a vortex line to move steadily together with superfluid. The kink supply rate 
required just before the glitches in the Vela pulsar is ~ 2 x 10 15 s -1 . This rate is much 
lower than the minimum frequency (equation (31)) of oscillations which can be excited in 
the pinned part of a vortex line. Hence, the oscillations of vortex lines cannot be a direct 
cause for the creation and supply of kinks. Moreover, energy of ~ 4.0 MeV (equation (23)), 
which is substantially larger than the thermal energy, is needed locally to form a kink. 
Energy source for kink formation is also pointed out as a difficulty included in the kink 
propagation model. We note that the kink solution (18) can also be expressed as a sum of 
fourier components of different wave numbers. As seen from the dispersion relation (30), the 
phase velocity of a vortex wave depends on wave number. The vortex waves with different 
wave numbers propagate with different phase velocities. Hence, even if a kink is formed 
around the end of a single crystal, the kink feature will be smeared out during propagation. 
Furthermore, our perturbation analyses show that the vortex configurations with static kink 
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structures are stable. These facts indicate that the kink motion following equation (32) is 
less likely. If the kinks do not move along a vortex line with the velocity given by equation 
(32), the conclusion by Jones is not guaranteed. 

We show in the previous sections that the vortex equilibrium configuration is in general 
composed of the static kink and straight segments. A vortex line in equilibrium lies deep 
in the pinning potential well and is strongly pinned to the lattice nuclei in its most part 
except for the kink part. From the perturbation analyses we find no unstable mode that 
grows with time for the pinned equilibrium configurations both with and without kinks. 
The pinned equilibrium configurations, though oscillations can be excited, may be stable. 
Hence, the static kink is a stable structure. Our studies show that strong pinning of vortex 
lines is likely to occur especially when the average direction of a vortex line is parallel or 
approximately parallel to the major axis of a crystal lattice and the Magnus force is not 
so strong. The orientation of a crystal lattice varies from one domain to another in a 
polycrystalline structure. The pinning strength decreases in the domain where the major 
axis of a crystal lattice is not close to the vortex line. However, as long as the equilibrium 
vortex configurations are represented by the kink solutions, pinning of a vortex line in a 
polycrystalline structure would be still strong enough to explain the large glitches (see the 
subsection 3.7). 



6. CONCLUDING REMARKS 

We have shown that the equilibrium configuration of a vortex line consists of straight 
and kink parts when the rotation axis of the star is approximately aligned with the major 
axis of a crystal lattice and the Magnus force is weak. The straight part is much longer than 
the kink part. Furthermore, the straight parts pass through nearly all of nuclei along the 
major axis, whereas the kink parts lie mostly in space between the lattice planes. The vortex 
line is strongly pinned to the lattice nuclei in the straight part, while not in the kink part. 
We have found from the perturbation analyses that there exist two modes, rotational and 
helical (Kelvin), for oscillations excited on a vortex line. The vortex oscillations are possible 
only above a minimum frequency (Link, Epstein, & Baym 1993) and the phase velocity 
of oscillations varies with wave number. We find no unstable mode that grows with time. 
These facts lead to the conclusion that the vortex configuration with kinks may be stable 
and yield the strongly pinned state especially when the vortex line is close to the main axis 
of a crystal lattice. 

Jones (1997, 1998) argues that pinning is not strong enough to explain the large glitches 
adopting a kink motion in a polycrystalline structure. Our studies show that the static 
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kink is a stable structure and the kink motion is less likely. If the kinks do not move 
along a vortex line, the conclusion made by Jones is not guaranteed. Nuclei in the crust 
may be in a polycrystalline structure. The orientation of a crystal lattice varies from one 
domain to another. The pinning strength decreases in the domain where the major axis 
of a crystal lattice is not close to the vortex line. However, as long as the equilibrium 
vortex configurations are represented by the kink solutions, pinning of a vortex line in a 
polycrystalline structure would be still strong enough to explain the large glitches. 

We have shown that the Magnus force makes the vortex line parabolic on global scale 
and causes it to incline against the major axis of a crystal lattice. As the Magnus force 
increases and the rotation axis inclines against the major axis of a crystal lattice, we need 
to include the major axis components for force and displacement, though neglected in the 
present study. We expect that the lattice planes different from those considered here may 
become important for pinning. 

The pinning sites are discrete although we have assumed the uniform pinning potential 
along the major axis of a crystal lattice. Link, Epstein and Baym (1993) find that there 
appears the band structure in the energy spectrum for excitations on a vortex line when 
the discrete nature is included in the pinning site distribution. It is also necessary to take 
account of the discrete nature of the pinning sites when we consider the vortex equilibrium 
configurations with the scale length comparable to the lattice constant. 

The above arguments suggest that the three dimensional calculations adopting the dis- 
crete pinning sites should be conducted in near future in order to understand the equilibrium 
vortex configurations and vortex oscillations and judge the relevance of vortex pinning for 
the pulsar glitches correctly. 

We thank M. Ruderman, R. Epstein and B. Link for many useful discussions. We also 
thank S. Nakamura for his help in calculations. This research was supported in part by the 
Grant-in Aid for Scientific Research (C) (10640234, 12640229, 12640302). 
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Fig. 1. — A straight line solution to a vortex equilibrium configuration. The major axis of 
a single crystal is chosen as a z-axis. The horizontal coordinate expresses displacement of a 
vortex segment from the z-axis in the direction of a reciprocal lattice vector. All coordinates 
are multiplied by g = 27i/a. This straight solution is obtained when the rotation axis of 
the neutron star is perfectly parallel to the z-axis. Here the Magnus force is included. The 
vortex line passes through almost the bottom of a pinning potential (0 = — ip ~ 0.016). 
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Fig. 2. — Sinusoid-like solutions to a vortex equilibrium configuration. The Magnus force 
is taken into account. The average orientation of a vortex line, which corresponds to the 
direction of a rotation axis of the star, is parallel to the z-axis. The integration constant in 
equations (16) is larger in the dashed line than in the solid line. 
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Fig. 3. — An equilibrium configuration of a vortex line with one kink, derived with no 
Magnus force. The vortex moves from one lattice plane (0 = 0) to the adjacent one (0 = 2tc) 
by forming a kink. Except for the kink part, the vortex line is straight and parallel to the 
z-axis. 
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Fig. 4. — Vortex equilibrium configurations with kinks for different inclination angles of a 
rotation axis against the z-axis. Here the Magnus force is ignored. The numbers on the 
curves denote the inclination angle. The average orientation of a vortex line corresponds 
to the direction of a rotation axis of the star. The main lattice planes are located at the 
horizontal coordinate of <fi = 2nir. As the rotation axis inclines against the z-axis, the pinned 
straight lengths decrease and kinks become less prominent. 
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Fig. 5. — Vortex equilibrium configurations when the Magnus force is included. The su- 
perfluid flows in the y-direction. The numbers attached to the curves denote the superfluid 
velocity relative to the lattice nuclei. The boundary conditions adopted are <fi = at z = 0, 
and d(j)/gdz = 1.4 x 10~ 3 , 6.6 x 1CT 3 , 2.4 x 1CT 2 at z = respectively for the cases of 
v s = 10 2 , 10 3 and 10 4 cm s -1 . The parabolic structure appears on the long scale and its 
curvature increases with the relative velocity. As the Magnus force increases, the pinned 
straight lengths decrease and the kinks become less prominent. 
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Fig. 6. — Dispersion relations for oscillations allowed on a vortex line. The numbers by 
the curves denote a displacement of the vortex segment, <f> , from the z-axis in equilibrium 
configuration. The oscillations in the pinned segments are possible only above a minimum 
frequency. 



